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Abstract. We demonstrate that a Langevin equation that describes the 
motion of a Brownian particle under non-equilibrium conditions can be exactly 
transformed to a special equation that explicitly exhibits the response of the 
velocity to a time dependent perturbation. This transformation is constructed 
on the basis of an operator formulation originally used in nonlinear perturbation 
theory for differential equations by extending it to stochastic analysis. We find 
that the obtained expression is useful for the calculation of fundamental quantities 
of the system, and that it provides a physical basis for the decomposition of the 
forces in the Langevin description into effective driving, dissipative, and random 
forces in a large-scale description. 
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1. Introduction 

Owing to technological advances in the methods of manipulating systems on sub- 
micrometer length scales in aqueous solution, there is increased interest in studying 
the nonequilibrium nature of such small systems. Recent studies have yielded several 
universal relations, including the fluctuation theorem U, the Jarzynski equality [2] 
and the Hatano-Sasa nonequilibrium steady state equality and the validity of 
these relations has been verified experimentally in small systems consisting of beads 
and bead-RNA complexes ^JEIEI- In addition to the verification of these equalities, 
there is an interesting experimental study providing new information concerning the 
nature of nonequilibrium systems, in which the velocity response to an external 
perturbation and the correlation of the corresponding fluctuations were measured in 
a bead system |7] . We should also mention that there has been a substantial progress 
in understanding stochastic resonance 

Although great progress has been made in the experimental investigation of such 
small systems, direct measurements have been possible only for kinematic quantities, 
such as the position of a bead. Given this situation, in order to elucidate the mechanics 
of small systems, it is necessary to extract information regarding mechanical quantities, 

§ To whom correspondence should be addressed (harada@chem.scphys.kyoto-u.ac.jp). 
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such as 'force' and 'energy', from kinematic quantities. However, this task is not easily 
accomphshed in general For instance, consider the problem of guessing the equation 
of motion of a motor protein on the basis of experimental results. It can then be 
understood why it is difficult to determine an effective potential for the center of mass 
and to evaluate the size of dissipation effects, because neither a canonical distribution 
nor a fluctuation-response relation exist for nonequilibrium systems. (By contrast, 
these can be utilized to determine an effective potential and a dissipation strength 
near equilibrium states.) 

In this paper, we study a nonequilibrium Langevin system for the purpose of 
developing a new theoretical method to extract information regarding mechanical 
quantities in small systems. Specifically, we attempt to re-express a Langevin equation 
in terms of directly measurable quantities. Here, let us recall that for differential 
equations, there exists a theoretical method by which, for example, a Rayleigh equation 
that describes nonlinear oscillations can be transformed perturbatively into a simple 
equation displaying the observed frequency ^Ul- Following this idea, it might be useful 
if we could transform a Langevin equation into a special form that makes mechanical 
properties manifest. Specifically, motivated by a recent phenomenological study of 
energy efficiency we regard the response to a time-dependent perturbative force 
as a key property that connects kinematic and mechanical quantities. 

Here, we point out that a Langevin equation is equivalent to the corresponding 
Fokker-Planck equation as long as we consider statistical properties. However, when 
we study problems related to force and energy, analysis of a Langevin equation is more 
appropriate, because it can be regarded as an equation of force balance. Furthermore, 
physical quantities for a single trajectory, which are often observed in bio-mechanical 
systems jl2) . can be described only by a Langevin equation. For these reasons, we seek 
a useful representation of a Langevin equation rather than a Fokker-Planck equation. 

Explicitly, we study the simple one-dimensional Langevin equation 

^i{t) = f-u'{x{t))+m+£Mt)- (1) 

Here, 7 is a friction constant, / is a constant driving force, U{x) is a periodic potential 
of period £, and the prime denotes differentiation. Also, ^(t) is Gaussian white noise 
satisfying 

mm)^^iTS{t-s), (2) 

where (• ■ •) denotes the average over all noise histories. (See the references ^2) and 
for discussion of the physical basis of this equation in the case e =0.) The term 
sfpit) with sufficiently small e represents a perturbation force we use to investigate 
the response of the system. An initial condition x(to) — xq is given at t — tg, and we 
often take the limit tg —^ —00. In the argument below, we denote the average over all 
noise histories in the limit to —00 as (• • ■)^. Under these conditions, there exists 
the time-dependent distribution function poo{0,t), such that the relation, 

{A{xm^^ f d0A{e)p^{0,t), (3) 
Jo 

is satisfied for any ^-periodic function A{9). Note that poo{0,t) can be expanded in 
the form 

Poo{0,t)=Pst{0)+ep'^'H9,t) + Oie^). (4) 
With the above preparation, in Section |21 we prove the following theorem. 
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Theorem : Under the condition to — oo, the equation can be transformed into 
the form 



dsL{s){x{t - s) - vs) = / dsL{s)at^s)+m+£fpit), (5) 
Jo 

where Vs = {x{t))g is the steady state velocity. The function L{s) here is determined by 

the system parameters 7, /, U{x) and T, and is independent of e. L{s) is determined 

so that it satisfies the causahty property, i.e., L{t) = for t < 0. The function ({t) 

depends on {x{s)}ta<s<t, {£.{s)}t„<s<t, and {/p(s)}to<s<t, and most importantly it 

satisfies the relation 

(CW)e = 0{S'). (6) 

These functions, L{s) and C{s), are determined from the set of eigenfunctions of the 
Fokker-Planck operator corresponding to the Langevin equation as explained in 
Section O 

We now explain the physical significance of this theorem. First, note that (0) is 
equivalent to the form 

/•oo 

x{t)-v, = at)+ dsR{smt-s)+efp{t-s)). (7) 







Here, the quantity R(t) is determined from L(t) through the relation 

i?H = Lito)-\ (8) 

where R{lo) and L{uj) are the Fourier transforms of R{t) and L{t). R{t) for t > is 
obtained by the inverse Fourier transform and we also require the causality as R{t) = 
for t < 0. Throughout this paper, we use similar notation, with the Fourier transform 
of an arbitrary function A(t) given by 

/oo 
dte'"*A(<). (9) 
-00 

Then, the average of Q) becomes 

POO 

{x{t))^ -v,^e I dsR{s)fp{t ~s) + 0{e^). (10) 



"'0 

This equation characterizes the linear response of the velocity to the perturbation. 
We therefore call R{s) a "response function". Its Fourier transform, R{uj), is called 
the "dynamic susceptibility". The expression ||7J), which is equivalent to the original 
equation of motion , explicitly represents the response to a perturbation. 

Next, we show how the expression (jSJ is related to the problem of force 
decomposition in a large-scale description. Let us express the force acting on the 
particle at time t as 

cj,{t) = f~U'{x{t)) (11) 

and consider a finite time average of (/'(t). For convenience, we express the finite time 
average of a quantity A{t) as 



1 



t + T 



Ar{t) = - / dsA{s). (12) 
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From |(TJ, ©, l(ni) and (jHI, we find 



^ ^^i^+lL^ + / dsL{s)Cr{t 











We now assume that r is sufficiently larger than the inverse of the smallest decay rate 
of L{s). Then, because Ar{t) varies slowly with respect to t, (|13|l becomes 

Mt) ^ mvs ~ (i(o) - ^)^(i±il^^ + mcrit)- (14) 

T 

This expression implies that a time averaged force '/>r(i) can be decomposed into a 
driving force (the first term), a dissipative force (the second term), and a random 
force (the third term). In the reference tl5J, the force decomposition for a large-scale 
description with Q is investigated, and it is shown that the condition 

lim t(C.C.)o = (15) 

uniquely determines the force decomposition of the type given in 1)14(1 . In other words, 
only the component orthogonal to is regarded as a random component of 0t- in such 
a description. In Section [2. 41 we confirm that this orthogonality condition is satisfied. 
Actually, we have arrived at the expression (jsj by seeking a definition of C,{t) for 
which 1(15(1 holds. In this sense, 1(15(1 was the guiding principle used in obtaining the 
transformation yielding Q from . 

In Section El we demonstrate that the expression is useful for calculating 
fundamental statistical quantities. As already stated, the susceptibility R{uj) is 
obtained as the Fourier transform of R{t). In [Appendix A| we provide a consistency 
check of our theory by comparing this quantity R{lli) with an expression that is 
obtained directly from the definition of the dynamic susceptibility. Furthermore, the 
Fourier transform C{lu) of the time correlation of velocity fluctuations, 

C{t)^{{x{t)~v,){x{0)-Vs))o, (16) 
can also be calculated directly from {Tj). In particular, for the case a; = 0, we can 
derive the following formulae: 

' (j^^-^-C)) 
' " (/>/-«-))' ' 



where 



I±{e) = / d(?'e±/3t/(e)T/3t^(fTe')-/3/e'^ (19) 





with P = 1/T. The key technical lemma used in the derivations of the above relations 
is proved in [Appendix B| 

We note that R{0) is equal to the differential mobility, dvs/df, and that the 
expression ((17(1 was calculated from the functional form of Vs{f) given in the reference 
[16[. Also, (7(0)72 is equal to the diffusion constant D defined by 

D^l\m^^{{x{t)^x{0)~vstf)^. (20) 
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The expression ()18f) was first derived in the reference ,14 employing a recursion 
formula that is established in the theory of stochastic processes. Subsequently, the 
same expression was obtained in a perturbation theory treatment applied to the 
Fokker-Planck equation corresponding to (Q) The method of derivation used 

in the present work is different from those used in the previous studies. 

Furthermore, from the correspondence of i?(0) with dvs/df and from the definition 
of D in l|2()(l . the nature of each term in the force decomposition H14|l is clearly 
understood. By substituting H14(l into the r-averaged form of we obtain a large- 
scale description of the Langevin equation: 

~^^^^ x{t + T)-xit) ^ ^^^^^^ ^ -^^^^ ^ mCrit)- (21) 
T 

The effective dissipation coefficient, L{Q), is the inverse of the differential mobility R{0) 
as explained in The magnitude of the fluctuating forces, ((fr(i) + ^(0)Ct(0)^)o 
is expressed in terms of the diffusion constant as 2D / R{0)^ , which provides the 
renormalization of the bare noise intensity 2jT. These quantities are obtained by 
use of the explicit expressions given in (ITTjl and ((TH|l . 



2. Proof of Theorem 



2. 1 . Operator formulation 



In order to make the argument mathematically clearer, we express the Langevin 
equation ^ as 



-idxit) = {f- U'{x(t))} dt + y^2Yrdw{t) + efpit)dt, (22) 

where w{t) represents a Wiener process JJ7j . Let A{x) be an arbitrary function of x. 
Then, Ito's formula gives 

dA{x{t)) = EM^A\x(t))dt + -A"{x{t))dt 

7 7 

rxr 1 

+ ^ / — A'{x{t)) ■ dw + -efJt)A'{x(t))dt, (23) 
V 7 7 

where 

F(x) = / - U'{x). (24) 

Here, • represents the product in the Ito interpretation. In conventional notation, this 
is written as 

'^^^ = AA(x)|_,(,) + A'{x{t)) ■ + £/p(t)), (25) 

where 

A=^MA + Z^_^ (26) 
J dx 7 dx"^ 

We can express the solution of (|25|l in the form 

A{x{t))^g{t)A{x)\^^^^^, (27) 

by introducing a time-dependent operator, g{t), that does not depend on A{x). 
Substituting H27I) into (|25() . we obtain an equation for Q{t) as 

^ = git)A + g{t)^ ■ + ef^it)). (28) 
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With the initial condition Q{tQ) - 



1 , the formal solution of this equation is derived as 
d 1 



dsQ{t — s)-^e 



ax 7 



s)+e/p(i-s)).(29) 



Thus, because the force acting on the particle at time t is given by 
= F{x{t)) 
= G{t)F{x)l^,^, 



(30) 



we obtain 



0(t) 



- Jt-tQ 



t-t„ 



X — Xq 

^ ^{t~s)+efp{t~s) 
7 



(31) 



x—x(t — s) 



where we have used H27|l in the derivation of the second term on the right-hand side. 

In this way, we have reformulated the original nonlinear stochastic differential 
equation in terms of the operator A. We note that a theoretical framework for 
carrying out such a reformulation in the case of differential equations was developed 
through application of nonlinear perturbation theory |18j . Also, analysis based on 
the "microscopic distribution function" |19| is essentially the same as the present 
formulation. 



2.2. Preparation of functional space 

In order to make our investigation of (|31|l concrete, we introduce a functional space, 
7i, consisting of all complex valued, square integrable, periodic functions of 6 on the 
interval < 9 < £. We endow this space with the inner product 

{h,,h2) = I dehi{e)h2{e), (32) 



Jo 

for /ii,/i2 £ Ti., where * denotes complex conjugation. All the eigenvalues, —Xj, and 
the corresponding eigenfunctions, ^j{d), of the operator A in H are determined by 
the equation 

A^,{9) = -X,<^>,i0), (33) 

where the index j — 0,±1,±2,--- is determined by the relations \j — Al^ and 
Ao = < Re(A±i) < Re(A±2) < • • • holding among their eigenvalues. When a 
complex eigenvalue happens to be degenerate, the corresponding labeling is modified 
appropriately. 

Because A given in H26II is not a Hermitian operator on 7i, it is convenient to 
introduce the adjoint operator of A through the relation 

(£fp/^i,/i2) = (/^i,A/i2), (34) 
where £fp is the Fokker-Planck operator: 

\ d T 

Note that the set of eigenvalues of £fp is identical to that of A. Then, we denote the 
eigenfunctions of £fp by (9) and choose their labeling so that we have 

CFP^,ie)^-xp,ie). (36) 
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We can choose these eigenfunctions such that the following hold: 

f de-^m^m ^5,,, (37) 

oo 

^ ^*{6)<i>,{e')^5{e-6'). (38) 

3 = -OQ 

In particular, we determine the normalization factor for the zero eigenfunctions 

*o(e)=ftt(e), (39) 
$o(e) - 1. (40) 
Because F{9) G we can expand F{9) in the form 

oo 

E (*.,^)'i>.w- (41) 

j = -oo 

Here, the following relation is easily confirmed: 

{■^o,F)=-fv,. (42) 
Then, from H31|l and (|41|l . we obtain 

0(0 = -tv, + bit) 

ft-to 

+ ds{at-s)+efp{t-s))-K{s,x{t-s)), (43) 

Jo 

with 

oc 

6(i)^ E (*.,J^)e-'^(*-*°^$,(xo) (44) 

and 

i^(s,0)^i E e-^^^(vI/,,F)$;.(0), (45) 

i=-oo 

for s > and i^(s, 6*) = for s < 0. Note that 6(t) satisfies the relation 

lim b{t) = 0. (46) 

^.5. Key step 

We decompose iir(s, 9) into two parts as 

K{s,9)^Ko{s)+K^{s,9), (47) 

where we have 

Ko{s) = {^o,K{s,-)), (48) 

for s > and Ko{s) = for s < 0. We choose this decomposition in order to satisfy 
the relation 



{K^{s,x{t)))^ = I d9pst{e)K^{s,9) 

(49) 



Fluctuating response equation from a Langevin equation 8 

Indeed, this equality can be derived from (|37|) and H39|) . Now, defining the quantity 

C{t) = - f " dsi^t - s) + efpit ~ s)) ■ K^{s, x{t - s)), (50) 
1 Jo 

we rewrite 0(f) as 

0(i) = -fvs + b{t) + f " dsKo{s){^x{t - s) - 0(t - s)) + 7C(t). (51) 



Here, we introduce 

^{t) =c^{t)--fVs, (52) 
A{t) EE x{t) - (53) 

and take the limit to —oo. The Fourier transform of (|51|l then yields 

^(w) = KoiLo){jA{Lo) ~ ^{lo)) + 7C(w). (54) 

From this, we derive 

^('^) = ^ i— [^o(^)7A(^) + (55) 

The inverse Fourier transform of this expression yields 

(j){t) ^ ^x{t) - dsL{s){x{t ~ s) ~ Vs) + dsL{s)C{t-s), (56) 
Jo Jo 

with 

i(t) ^ r l^^e— * 2 (57) 

for f > and L{t) = for f < 0. Then, recalhng (^TJ, we find that (jSEJ leads to the 
expression 

Next, using Q and we find that (^{t) defined by H5U|) satisfies 



(C(f)), = - / ds/p(f - s) {K^{s, xit - s)))„ + 0(£^), (58) 
7 Jo 

where we have used the Ito interpretation. Finally, from H49() . we derive ©. 
Remark 

As seen in the above discussion, the proper definition of Q is essential to obtain @. 
We arrived at the definition H5U|) in the following way. Intuitively, C, corresponds to 
what is in some sense the "random" part of the force 0(f) for the case £ = 0. With 
this in mind, we assume the following conditions: 

(C(i))o = 0, (59) 

{mm), = 0. (60) 

The condition H59|) is obviously necessary, and the condition H6U|I is inspired by the 
orthogonality condition H15|l reported in the reference ^^l- From the definition H50() . 
it is readily confirmed 

{amt'))o = - / ds {K^{s, x{t - s))), {at - s)m) 

I Jo 

= 0, (61) 

where H49f) has been used in the last line. It should be noted that the condition (|60|l 
is satisfied regardless of the interpretation of the multiplication between ({t) and ^(f). 
Then, guided by the condition (|60|l . we find the definition l|5UI) . and this leads to the 
proof of 
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3. Application 



First, we consider the response function R{t). From H45|l . 148|l and (|57|l . we obtain 
the following expression of the dynamic susceptibility: 



1 

1 + - 

7 



Xj — iuj 



(62) 



Replacement of the infinite sum in the right-hand side of (|62|l with a finite sum enables 
us to calculate R{uj) within a required precision because the contribution of the term 
with large \ j\ to the right-hand side in (|62|l becomes on the order of for sufficiently 
large The response function R{t) is obtained from the inverse Fourier transform. 

Next, note that the most important quantity characterizing the statistical 
properties of fluctuations described by Q might be the velocity correlation function 
H16|l . For the equilibrium case, / = 0, as a result of the detailed balance condition, 
C{t) is determined by R{t) through the fluctuation-response relation. However, for 
nonequilibrium cases, there is no such relation. Thus, we need to derive an expression 
of C{t) independently of R[t). 

Substituting 15U|) into {Tj), we find that it is convenient to define 



M(s,e) 



1 

-K^{s, 
7 



R{s) 



We then obtain 



C{t) = 27T / ds {M{s, e)M{s - t, 6l))o , 
Jo 

by use of Noting that C{t) = C{-t) and that M{t,e) 
C{u)) = {\M{uj,e)\^ 
( 

^-0 



2r 

7 



V 



1 



A, 



(63) 
(64) 

= for t < 0, we calculate 
(65) 



With this expression, Ciio) can be calculated with a required precision for the same 
reason as in the case of R{uj). The correlation function C{t) is obtained from the 
inverse Fourier transform. We remark that the result 165|) does not depend on the 
manner in which we interpret the multiplication between x{t) and i(0) in l|16|l . 

Here, we derive more compact expressions for i?(0) and (7(0). First, for later 
convenience, we define 



(66) 



Vi9) = U{9) - fe. 
Obviously, F — ~V' . Then, it is easy to confirm the equality 

e^f' /''d0'e-'3t'(»-e')-««' 
Jo 

= f d0'e-^^(«') - (1 - e^/^ df?'e-^^(^') . (67) 
Jo Jo 

Using this, the expression of Pst{S) obtained by solving the equation C-pppst{d) — 

with standard techniques (see, e.g. [21) can be transformed into the form 



Pst{e) = 



(68) 
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The equahty (|67|l also allows us to prove the relation 

(The proof of (|69|l is given in [Appendix B| ) Then, substituting (j68|l and into 
both linSi and with w = 0, we obtain (UTj) and ifT^ . 

Before ending this section, we discuss the similarities and the differences between 
our formulation and the standard analysis on the basis of the Fokker-Planck equation. 
We first emphasize that calculation of statistical quantities by use of (0 is in principle 
equivalent to that using the Fokker-Planck equation that corresponds to the Langevin 
equation In practice, we should use either formulation depending on the character 
of the concerned quantity. For example, statistics of an arbitrary function of x{t) can 
be more easily calculated in the Fokker-Planck formulation. A class of quantities 
such as the multiple-time correlation functions of x{t) — x{0), which is calculated in 
the Fokker-Planck formulation, can also be facilely calculated from by considering 
products of the x variables and subsequently averaging over statistical realizations. 

On the other hand, as an advantage of our formulation, the calculation of 
statistical quantities involving x{t) is facilitated by use of our formulation. For 
example, as seen in the above argument, the time correlation function of the velocity 
fluctuations, C(t), can be easily calculated in our formulation, while such a calculation 
seems rather difficult, but not impossible, in the Fokker-Planck formulation. Also, 
owing to the form of O, it provides a powerful tool to find a statistical quantity 
connected to the response function. Recently, an illuminating example using this 
advantage has been demonstrated in the reference |20| . in which the steady heat flux 
has been expressed in terms of the violation of fluctuation-dissipation relation. 



4. Concluding remarks 

In this paper, we presented the exact transformation of the Langevin equation to 
(jSJ as an example. Our theory may provide a novel way of understanding the "force" in 
a Langevin system as well as a new method to calculate R{uj) and C{uj) simultaneously. 
Due to the general nature of the argument used in its derivation, our theory can be 
applied to various stochastic systems. For example, it is straightforward to study 
a Biittiker model |21| using our theory. For other models, with time-dependent 
potentials [221 12S1 l2Hj , some additional techniques must be designed in order to 
construct a special expression that explicitly exhibits a response function. We will 
report such studies in a separate paper. 

The connection of the theory presented here to experimental studies is most 
important. As an example, suppose that we have an experimental system consisting 
of a small particle exhibiting a fluctuating movement on a one-dimensional track and 
that the mechanics of the system is under investigation. Then, the function L(s) can 
be determined from the response to a time-dependent perturbation by using JH)) and 
pOfl . and the statistical properties of C{t) can be determined from the time correlation 
of x{t) by using (|Sn|l . and (|^ . Given this L{t) and the statistics of ({t), we 
can write an equation of the form (O where C(i) is replaced with artificial random 
noise satisfying the observed statistics. Such an equation can be considered as the 
equation of motion for this system provided that the history dependence of ({t) is 
negligible. Although, strictly speaking, C depends on the history of x, this description 
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may provide a good starting point for the construction of a phenomenological theory 

El 

The idea of re-expressing differential equations reminds us of the normal form 
theory, whose goal is to transform the equation in question into as "simple" a form as 
possible. The normal form theory was extended to stochastic processes, with the main 
focus on bifurcation problems |26l 1271 1^ . Furthermore, mathematical techniques for 
random normal forms have been developed |29j . However, as far as we are aware, the 
transformation of to lO has never been presented. 

In another related work, a generalized Langevin equation with a memory function 
was derived from a nonlinear Langevin equation by employing a projection operator 
method jJUl- In this method, functions of dynamical variables are projected onto 
a subspace of interest, by using the so-called Mori identity PO]- Then, the choice 
of the projection operator represents the essence of the problem when considering 
whether the obtained result is physically meaningful, as pointed out in the reference 
|31|. However, as far as we understand, there was no such discussion in the reference 
| 19j . Of course, whatever method we use, the important point is whether the obtained 
result is related to quantities and relations that can be measured experimentally. The 
determination of dissipative and other forces should be done in such a way that it is 
clear how they can be measured experimentally. 
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Appendix A. Verifying the formula H57|l 

Here, we confirm the validity of the formula (|57|l by directly calculating the response 
function from its definition l|10|) . To this end, we utilize a perturbation force, fp{t), 
that takes the form of a step function: fp{t) = for i < and fp{t) = 1 for t > 0. 
Then, we have 

{x{t))^~v^^e [ dsR{s) +0{e^). (A.l) 



Using lO, the left-hand side of this equation can be expressed as 

(i(t)),-«s= f dep^{9,t)^^^^^^^^ -V,. (A.2) 



The function Poo{0^t) satisfies the equation 

j-Poc{0,t) = (Up - If (t) — )p^i9,t), (A.3) 
ot -f oO 

with the initial condition Poo(^5 0) = Pst{&)- 

Substituting the expanded form Q into HA.3|I and extracting terms proportional 
to e, we obtain 

§-/'Ho,t) = Cppp^'\e,t) - ^/p(tKt(^)- (A.4) 
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Next, we write p^^\9,t) in terms of the eigenfunctions of Cpp: 

oo 

P^'\9,t)= (A.5) 

j = -oo 

Substituting this form into IIA.4|) . wc obtain 

pf\t) = -X,pf\t) + ^/pW(*o, (A.6) 
Integration of ljA.6|) with the initial condition p^^ [0) = yields 

PrW = -/pW(*o,*;) fdse-^^^. (A.7) 
7 Jo 



Using the results (jA.5|l and (|A.7p . we express (|A.2|I as 



e 



A, 



7 Jo 



+ 0(e2). (A.8) 
Comparing this with (jA.l|l . we find 

jR{t)^6{t)+ E -(*„i^)e-^^*(*o,$;). (A.9) 

j = -oo 

The Fourier transformation of this relation gives 

7i?M = l+ E -^;^^(*o, (A.IO) 
which is found to be equivalent to 
Appendix B. Proof of 

Differentiating with respect to 6, we obtain 

A*$; = -A,$;, (B.l) 

with 

Because ff^- is an eigenfunction of A^ , below we employ the notation 

$] = (B.3) 
The operator adjoint to A^, Cpp, is defined through the following: 

(4p/ii,/j2) = (/^i,A*/i2). (B.4) 
The explicit form of Cpp is 



FP 

1 , d T d'^ 
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The set of eigenvalues, —Xj, of this operator Cpp is identical to that of A*. The 
corresponding eigenfunctions are labeled such that they satisfy the equation 

^Ip^-] = -A***. (B.6) 
These eigenfunctions can be chosen so as to satisfy the orthogonality condition 

(B.7) 

Then, differentiating (|B.6|) with respect to 9, we obtain 



/:fp^']' - -A**]'. (B.8) 
Thus, applying (|B.7|I . we have 

*f - -'i'j (B.9) 

for the case j ^ 0, and = 1. 

Now, we define the quantity Q{&) as 

E ^^^f^W- (B.IO) 

Using HB.Sp . we calculate 

b}Q = -V". (B.ll) 

7 

Solving this equation, we obtain 

Q{B) = -1 + gie'3^('') /' d^'e-'^^C^') + gae^^^'l (B.12) 



The two constants q\ and (72 here are determined by the conditions 

Q(0) = Q(^), (B.13) 
(4'S,Q) = 0. (B.14) 



The condition (|B.13|) leads to 
1 

92 q\ 



3/3/^ - 1 7o 



d0'e-^^(«'). (B.15) 



Substituting this into l|B.12|) . we obtain 



1 + Q{Q) = qi 



cf'f' - 1 
9i 



(B.16) 

where we have used the identity ()67|l . Using l)B.16|l in (|B.14|I leads to the final form 
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